Direct view on the ultrafast carrier dynamics in graphene 
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^H Abstract 

> 

'^ The ultrafast dynamics of excited carriers in graphene is closely linked to the Dirac spectrum 

}~^ and plays a central role for many electronic and optoelectronic applications. Harvesting energy 

■^ from excited electron-hole pairs, for instance, is only possible if these pairs can be separated before 

O 

^^ they lose energy to vibrations, merely heating the lattice. While the hot carrier dynamics in 

^ graphene could so far only be accessed indirectly, we here present a direct time-resolved view on 

/^ the Dirac cone by angle-resolved photoemission (ARPES). This allows us to show the quasi-instant 

thermalisation of the electron gas to a temperature of more than 2000 K; to determine the time- 
resolved carrier density; to disentangle the subsequent decay into excitations of optical phonons 
and acoustic phonons (directly and via supercollisions); and to show how the presence of the hot 
carrier distribution affects the lifetime of the states far below the Fermi energy. 



Exploiting the unique properties of graphene in electronic or optoelectronic devices in- 
evitably involves the generation of hot carriers, i.e. electrons or holes with high energies 
compared to the Fermi energy [TH5]. Such hot carriers are thought to thermalise via the 
electron-electron interaction on a time-scale of ~ 30 fs |1H6], leaving the electronic system 
at a well-defined temperature that can be substantially above the lattice temperature. If 
an initial electron-hole pair had been created by photo absorption, the thermalisation can 
be accompanied by carrier multiplication, i.e. the generation of more than one electron-hole 
pair per absorbed photon [TJ |H] . The hot electronic system can then efficiently lose energy 
by the fast emission of optical phonons, but this channel is blocked for electrons with too 
low energies [^ 200 meV) or when the temperature of the optical phonons reaches that of 
the electrons P4T2] . Then a slow decay via acoustic phonons sets in. Indeed, this decay is so 
slow and inefficient that three-body collisions with defects and acoustic phonons, so-called 
supercoUisions, can become the dominant contribution rather than the exception [T314T5] . 
This hot carrier dynamics, summarised in Fig. [T^, has been studied intensely but so far 
the ultrafast decay could only be probed indirectly, for example through the change of the 
transmissivity in a pump-probe experiment. 

Time and angle-resolved photoemission spectroscopy (TR-ARPES) offers an unparalleled 
access to the detailed carrier dynamics in the spectral function of solids and has been suc- 
cessfully applied to study many systems (see for example Refs. |16l [T7]). It has the potential 
to directly provide the out-of-equilibrium single-particle spectral function, the statistical dis- 
tribution of the carriers and its time evolution. However, for graphene the application of this 
technique is hampered by the position of the Dirac cone at the K point, far from the centre 
of the Brillouin zone. In order to collect TR-ARPES data from the Dirac cone, the final 
state electrons have to have a sufficiently high fcy and this can only be achieved for photon 
energies higher than ^ 16 eV, significantly higher than commonly available in conventional 
ultrafast laser sources [T71 dH]. In order to overcome this hurdle, we have used the coherent 
high harmonics generated by laser light at the Artemis facility at the Central Laser Facility 
/ Rutherford Appleton Laboratory. Hot electrons in graphene are generated by a 30 fs pulse 
with a photon energy of hu =0.95 eV. This pump pulse leads to a transition from the valence 
band to the conduction band of graphene, as illustrated in Fig. [T^. The spectral function 
around the Dirac cone is then measured by another laser pulse oi hu = 33.2 eV, following 
the first pulse by a variable time delay. 



In order to closely approach the electronic properties of pristine, free-standing graphene, 
we use a sample of hydrogen-intercalated so-called quasi free-standing monolayer graphene 
(QFMLG) on SiC [T^ I2U] . The sample is slightly hole doped with a carrier concentration 
of 5 X 10^^ cm~^, placing the Dirac point 240 meV above the Fermi energy Ep. QFMLG is 
characterised by high crystalline quality, sharp photoemission lines, weak electron-phonon 
coupling [21] and an efficient decoupling from the substrate, such that even subtle many- 
body effects can be observed [22]. ARPES data from the Dirac cone of QFMLG are shown 
in Fig. [T]b for a negative time delay, i.e. before the excitation with the pump pulse. As a 
comparison, high-resolution ARPES measurements from the same sample are shown in Fig. 

Fig. Wd-i show TR-ARPES spectra taken at different delays after the pump pulse (see 
also Supplementary Movie SI) and Fig. flk-i the corresponding difference between these 
spectra and the spectrum before the arrival of the pump pulse. A clear increase (decrease) 
of the spectral intensity in the conduction band (valence band) is already discernible in the 
raw data and even clearer in the corresponding difference. Note that the strong intensity 
asymmetry between the two branches of the Dirac cone in every spectrum is a well-known 
interference phenomenon in ARPES from graphene [23] . 

Immediately after the pump pulse, one should observe a depletion (increase) of spectral 
intensity around the energy of the initial (final) states of the excitation, or at least a separate 
distribution for electrons and holes [H] but it has been shown that the electron-electron 
interaction leads to a thermalisation of the carriers on a time scale of ~ 30 fs, faster than 
the experimental time resolution [l]-[6] . If this is so, we should only ever observe an electron 
gas with a thermal Fermi-Dirac (FD) distribution. 

We have direct access to the distribution of the carriers using the procedure introduced in 
Fig. [2] For each delay time, we slice the image into momentum distribution curves (MDCs), 
i.e. the photoemission intensity at a given binding energy as a function of fc||. A stack of 
such momentum distribution curves for a time delay of 100 fs can be seen in Fig. |2^. Each 
MDC is then fitted with a Lorentzian peak, giving the peak position and the integrated 
intensity of the peak as a function of binding energy. The MDC intensity as a function of 
energy and time is shown in Fig. |2]b. It can be viewed as the statistical distribution of 
the carriers, tracking the dispersion of the state. The analysis of the MDC peak position 
confirms that there is no time-dependent change of the band structure, as expected for the 
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FIG. 1: a Dirac spectrum of graphene with an optical excitation hvpp and a probe hupB- Fol- 
lowing the excitation, electrons (black spheres) thernialise via impact scattering (magenta curved 
arrows) or Auger recombination (black curved arrows). Decay into holes (orange spheres) in the 
valence band can be mediated by emission of optical phonons (red wiggled line) with energy hvop 
or supercollisions involving acoustic phonons (green wiggled line) with energy hvsc recoiling on 
impurities (dashed green line), b Photoemission intensity around the Dirac cone for a negative 
time delay (before the pump pulse), c Equilibrium high-resolution spectrum from the same sample 
taken with synchrotron radiation. Insert: Brillouin zone with measurement direction (dashed red 
line), d -f Spectra taken at different time delays, g - i Corresponding difference spectra, i.e. change 
with respect to the spectrum before the pump pulse. 

low fluence employed here. Note that the apparent upward bending of the dispersion near 
Ep is most likely caused by the finite energy resolution |24j . 

The direct access to the electron distribution confirms the ultrafast thermalisation of the 
carriers. The data in Fig. [2)3 are well described by a FD distribution for all time delays, as 
illustrated by a few cuts in Fig. |2)i-k. This confirms that the thermalisation of the carriers 
happens at a time scale that is shorter than the experimental time resolution, permitting 
the direct determination of an electronic temperature Te(t). Interestingly, deviations from 
the FD behaviour can be observed immediately after the main pump excitation at time 
delays between and 200 fs where the photoemission intensity far above Ep is higher than 
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FIG. 2: a Spectrum from Fig. [T]d shown as momentum distribution cuts (only a sub-set is shown) 
and fits of the individual cuts by Lorentzians to determine intensity and position, b Intensity 
obtained in this way plotted as a function of energy and time delay, c Fitted MDC peak positions 
at various time delays, d - j Cuts through the data in b with fits by Fermi-Dirac distributions 
convoluted with a Gaussian to account for the energy resolution, k Fitted distributions at selected 
time delays. 



expected for the FD distribution. While this is a strong indication of an out-of-equilibrium 
situation before thermalisation, the limited time resolution does not permit its more detailed 
investigation. 

The electronic temperature extracted from this procedure is shown in Fig. [3^. The most 
simple analysis shows that the data cannot be fitted by one single exponential decay but 
by a double exponential decay involving two time constants, ti ^ 160 fs and T2 ~ 2700 fs. 
These are tentatively assigned to processes involving optical phonons and acoustic phonons / 
super collisions, respectively. In order to analyse this further, we employ a phenomenological 
three temperature model in a similar form as recently applied to a TR-ARPES study of 
a high-Tc7 superconductor [121 I2S]- The model assumes that three different temperatures 
exist in the system: the electronic temperature Te, the temperature of the optical phonons 
Tp and that of the acoustic phonons Ti. This leads to three coupled differential equations: 
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where Cg, Cp and Ci are the heat capacities, Ai and A2 are the electron-phonon couphng 
constants to the optical and the acoustic phonons, respectively, VtEm is the energy of the 
optical phonons, modelled as Einstein oscillators, rie and Up are distribution functions given 
as {(P-Einl^BT{^ei-p) _ i)-i^ g(fi(Te)) is the density of states evaluated at the temperature 
dependent chemical potential fi(Te), I is the mean free path of the electrons, Vs is the sound 
velocity in graphene and kp is the Fermi wave vector with respect to the K point. 

Details of the model are given in the Supplementary Material. Here we concentrate on 
the general form of the equations. The first term in ([I| stems from the heating by the pump 
pulse S{t) that has a Gaussian distribution with a full width at half maximum (FWHM) of 
30 fs. /3 is a constant with the dimension of a heat capacity but note that the heat capacity 
of the electronic system is not properly defined during the pump phase. After the excitation, 
the thermalisation of the electrons is assumed to be sufficiently fast that a definition of Tg is 
meaningful and Cg can be defined. The electrons can now lose energy through either coupling 
to optical phonons (second term in ([I])) or to acoustic phonons, either via supercoUisions [13] 
or directly [26] (remaining two terms). The corresponding equations for Tp and Ti describe 
how these sub-systems are heated up by the energy that is lost in the electronic system. In 
addition to this, it is also possible for optical phonons to undergo anharmonic decay into 
acoustic phonons with a characteristic lifetime Tp. 

The unknown parameters in the equation system are the coupling constants Ai_2 and the 
scaling parameter /3. By adjusting these, it is possible to obtain an excellent fit to the data 
(black line in Fig. pa). Note that (3 only determines the initial temperature, it has no 
infiuence on the curve after the pump pulse duration. We find that the coupling constants 
are very small with Ai = 0.052 and A2 = 0.0042, as expected for lightly doped graphene 
[271 [28] • Indeed, our value for Ai is consistent with the electron-phonon coupling strength 



determined in Ref. [12] and so is the fact that most of the incoming energy is absorbed by 
the electronic system. 

We can also test if it is possible to fit the curve in the absence of supercollision processes 
by removing the corresponding term. This results in the grey dashed line. It clearly shows 
that the supercollisions are essential once Tp = Te, the point at which cooling by optical 
phonons is exhausted. We also observe that the lattice temperature increases very little, 
which can be explained by its much higher heat capacity, allowing it to act like a perfect 
heat bath. 

Given the fact that the system follows a FD distribution with a well-defined temperature, 
we can calculate the time- dependent photo-induced carrier density rij, i.e. the number of 
photo-excited electrons above Ep (see Fig. [sja and details in the Supplementary Material). 
From rij and a calculation of the absorbed energy at our maximum electronic temperature 
we arrive at an estimate of the number of electron-hole pairs n' generated by the pump pulse. 
This allows us to quantify the carrier multiplication {CM) in the system, i.e. the number 
of electron hole pairs generated for every photon absorbed. In our case, CM remains below 
unity, merely reaching ^0.5, consistent with theoretical predictions [29]. Note that these 
predictions also show that CM can reach values well above unity for the lower fluence values 
and higher photon energies actually relevant for photoelectric applications. Here, because 
of the need for a minimum fluence for a pump-probe experiment, the CM > 1 regime is 
outside our experimental conditions. 

The linewidth of states in ARPES can give detailed information about the many-body 
effects in the hole spectral function. In particular, the MDC width can be related to the 
imaginary part of the self-energy or, equivalently, to the lifetime of the photohole. We can 
therefore not only investigate the ultrafast redistribution of carriers but also the effect this 
has on the energy and momentum-resolved lifetime of a given state. Fig. |4^ shows a cut 
through the spectral function deep in the occupied states, at a binding energy of 850 meV, 
for different time delays. After the pump pulse, not only does the total intensity decrease 
due to the lower value of the FD distribution, we also observe an increase of the MDC width 
from ^ 0.115 A^^ to ~ 0.141 A^^. This change in width corresponds to the opening of an 
additional decay channel for the photohole. We estimate the lifetime from the MDC width 
and find that the initial lifetime tq ^ 1.26 fs is reduced to ti ^ 0.92 fs. The characteristic 
lifetime for the new decay channel r^ can be estimated from 1/ri = I/tq + l/r^ to be r^^ ^ 
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FIG. 3: a Electronic temperature Te{t) as a function of time delay (markers). Fitting the data to 
the three temperature model provides an electronic temperature fit with (black line) and without 
(dashed grey line) the supercollision (SC) term. From the fit we also obtain the temperatures of 
the optical phonons (blue) and the acoustic phonons (red), b Photo-induced carrier density n/ 
and carrier multiplication calculated from the fitted electronic temperature. 



3.41 fs. Fig. |4]b shows the time-dependence of the photohole hfetime estimated from the 
MDC width, showing that the new decay channel of the photohole closes very quickly on 
a timescale of 200 fs. This observation permits us to narrow down the origin of the decay 
channel. Consider the possible hole- hole scattering processes shown in Fig. |4|3 [30] . Process 

1 involves scattering with holes at a lower binding energy than the photohole while process 

2 requires holes with a higher binding energy. Since the number of these high energy holes 
decreases very rapidly [31], the time-dependence of the MDC broadening suggests that the 
new decay channel is dominated by process 2-like scattering events. Note that the decay in 
Fig. |4^ appears very similar to the initial decay of Tg, a process ascribed to the excitation 
of optical phonons. However, optical phonons are not directly responsible for the additional 
decay channel of the photohole because such processes are always possible; they are not 
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FIG. 4: a MDC cuts through the spectral function at a binding energy of 850 nieV for different time 
delays, b Lifetime of the hole state as a function of time. The line is a fit to an exponential model 
with a time constant of 200 fs. Inserts: Dirac spectrum after excitation with hot electrons above 
and excess holes below the doping level. This situation gives rise to two hole-hole recombination 
processes (arrows). The dashed line marks the energy of the photohole. 



enhanced by the presence of hot holes. The reason for the time scale being similar is that 
process 2 requires high energy holes and these are rapidly depleted during the initial cooling. 
Concluding, we have been able to directly measure the time, energy and momentum- 
resolved statistical distribution of hot electrons in quasi free-standing graphene after a photo- 
excitation process. We confirm the ultrafast redistribution of carriers to a hot thermalised 
Fermi-Dirac function with no separate electron and hole distributions discernible within 
our time resolution. The decay process of this hot electron population proceeds with the 
expected involvement of optical and acoustic phonon excitations and we are directly able 
to reveal the role of supercollisions in this process. We quantify the number of induced 
electron-hole pairs and find that it stays under the threshold of carrier multiplication. This 



result is found to be consistent with quantitative predictions of carrier multiplication in 
graphene and suggests that this effect is likely to play a role in the low-fiuence situation of 
actual devices. Finally, we show that TR-ARPES can disentangle the photohole lifetime, 
described by the single-particle spectral function, from the lifetime of the excited charge 
carrier population distribution. 

METHODS 

TR-ARPES experiments were performed at the Artemis facility at the Central Laser 
Facility / Rutherford Appleton Laboratory [32]. A 1 kHz Ti:sapphire amplified laser system 
provided ultrafast (30 fs FWHM) infrared pulses at 785 nm, with an energy per pulse of 
10 mJ. For angle-resolved photoemission, 20 % of the laser energy was used to generate 
high-order harmonic femtosecond XUV pulses in a pulsed jet of argon gas. The time- 
preserving monochromator f33j selected the 21st harmonic with a photon energy of 33.2 eV. 
The remaining power was utilised to drive an optical parametric amplifier (HE-Topas), which 
provided tuneable laser pump pulses (30 fs). We chose a pump energy of 0.95 eV (1300 nm). 
The pump fluence was approximately 346 /xJ/cm^, and the beam was polarised perpendicular 
to the scattering plane (s-polarised) in order to avoid laser-assisted photoemission effects. 
For further details of the experimental set-up see Refs. [221 El]- Time, energy and angular 
resolution were set to 60 fs, 350 meV and 0.3°, respectively. The high-resolution static 
ARPES measurements in Fig. ^ were measured on the SGM-3 beamline of the synchrotron 
radiation source ASTRID in Aarhus |35j . 

H-intercalated epitaxial graphene on SiC(OOOl) was prepared ex-situ by the methodology 
given in Ref. [2H], and was cleaned by annealing to 550 K in ultra-high vacuum in order to 
remove absorbed water. The sample was held at room temperature throughout the entire 
experiment. 

We note that the extraction of the carrier distribution used in Fig. |2]is somewhat cumber- 
some compared to the frequently used method of analysing a cut through the photoemission 
intensity at a single angle or fcy. However, in the present case this simpler method leads to 
very inaccurate results for the electronic temperature, underestimating it by up to 1000 K. 

For the three temperature model in equations (|l|-(^, the following parameters were used: 
The optical phonons were approximated by Einstein oscillators with QEin = 200 meV, and 



10 



the anharinonic decay was set at r^ = 2.5 ps (see Ref. [H]). The mean free path of the 
electrons / was determined by high-resolution ARPES measurements from the MDC full 
width at half maximum at Ep, which is given by 1//. 

The lifetimes r in connection with Fig. |4] were also determined from the MDCs width 
and the slope of the band v = dE{k)/dk, using that / = vt. In the actual determination of 
r, however, the finite momentum and energy resolution have also been taken into account. 
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